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ARRANGEMENT OF BODIES OF REVOLUTION IN SUPERSONIC 
FLOW TO REDUCE WAVE DRAG 
By Morris D. Friedman 


Summary 

The wave drag of a combination of slender bodies of revolution at 
zero angle of attack is studied with a view to determining the arrange - 
ments for which the total drag is a minimum. Linearized theory is used 
to calculate the pressure distribution in the field surrounding the 
bodies. The interference drag coefficient is computed for different 
arrangements . 

The special cases of two bodies and of a three-body combination with 
bilateral symmetry are considered. The bodies treated are of the form 
determined by Sears and Haack as having minimum wave drag for prescribed 
volume and length. They also have equal fineness ratios. Numerical 
calculations of the drag coefficient of interference are carried out and 
curves are drawn which show the relative positions at which mi nimum drag 
occurs. 

A three-body configuration is found for which the total wave drag 
is about 35 percent less than the sum of the individual wave drags of 
the three bodies. 

INTRODUCTION * 


The drag of a body of revolution in supersonic flow has been con- 
sidered by von Karman, Haack (references 1 and 2), and others. Haack 
and Sears (references 2 and 3) have determined theoretical body shapes 
for which the wave or pressure drag is a minimum. Such bodies have 
important present-day flight applications since this wave drag, as an 
additional form of drag at supersonic speeds, limits the performance of 
modern aircraft. 

If, in respect to the wave drag, the components of an aircraft such 
as the fuselage and wing-tip tanks and nacelles, say, are replaced by 
equivalent bodies of revolution, then the question of the combined wave 
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drag is of interest. Different arrangements of the bodies of revolution 
at zero angle of attack are investigated and that combination which has 
the smallest combined wave drag is determined. A practical combination 
would consist of a large body or fuselage and/or two smaller bodies 
which could represent either wing nacelles or tip tanks. In particular, 
the drag coefficient for the case of two slender bodies of equal fine- 
ness ratio (ratio of total length to maximum diameter) and unequal 
lengths is considered. (See fig. 1.) The procedure to be outlined, 
however, permits generalization, for the body shape chosen, to any num- 
ber of bodies of any relative size. 


SYMBOLS 


a speed of sound 

b abscissa of center of body of revolution 

c constant, the value of which determines the fineness ratio 

Cdj coefficient of interference drag 

Ci>j. drag coefficient of the configuration based on the total fron- 

tal area of the configuration 

Cp pressure coefficient 

Di drag change due to interaction of body pressure fields 

E complete elliptic integral of the second kind of modulus k 



F ° £ k 

K complete elliptic integral of the first kind of modulus k 

k modulus of the complete elliptic integrals 

L path of integration 

m half-length of body 

M Mach number 




RESTRICTED 


! 


NACA RM A51I20 


RESTRICTED 


<lo free-stream dynamic pressure 

r cylindrical coordinate (J y 2 + z 2) 

r Q particular value of r 

R • radius of a section 

Rttiay radius of maximum section 

S cross-sectional area of body or bodies 

V free -stream velocity 

x,y,z rectangular coordinates 

B • JW - 1 

C argument of the elliptic integrals of the third kind 

Ao term which occurs in the circular case of the elliptic integral 

of the third kind 

£ dummy variable of integration 

<P perturbation velocity potential 


Subscripts 


1 parent body 

2 satellite body 

I, II, III regions of integration 
t combination as a whole 


GENERAL CONSIDERATIONS 


The type of configuration studied is illustrated in figure 1. It 
consists of a large and a small body each of the form having minimum 
wave drag for given volume and length and each of the same fineness ratio. 
The large body is situated on the axis with its center at the origin. 

The smaller body has its axis parallel to the axis of the large body and 
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may be. within the flow field of the large body. The equation for the 
shape of such bodies, which may be supposed to be created by distributing 
sources and sinks along an axis, is given by Haack (reference 2) as 


R 2 = c [m 2 - (x-b) 2 ] 3/ ^ 2 


( 1 ) 


with b = 0 for the parent and b = b for the satellite body. (For 
b = 0, this formula is the expression for a thin body of revolution of 
length 2m with center at the origin and fineness ratio 2m/2R max .) 

To calculate the interference drag of such a configuration, two 
possibilities must be considered: 

1. Only one of the bodies is within the disturbed flow field of 
the other. 

2. Each body is in the disturbed flow field from the other. 

\ x 

When one of the bodies lies within the disturbed flow field of the 
other, the effects of the flow field in which this body is located must 
be considered. In other words, if all or part of either body is behind 
the Mach wave from the nose of the other body the streamlines will be 
distorted and a pressure will be exerted by the flow field of, one body- 
on the other. •••••' ' - 

The potential field which results from the interaction of the flow 
fields consists of the sum of the individual potentials and an interfer- 
ence potential. In the cases when there are multiple reflections 
between the bodies a series of interference potentials may occur. Since 
the interference potentials are usually of higher order of smallness, 
they will be assumed negligible. Of course, at very high Mach numbers, 
or when the bodies are relatively close to each other, the effect of 
this interference may not be negligible. 

The interference drag coefficient due to the location of a body in, 
the flow field of another can be evaluated by integrating the product of 
the additional, disturbed, pressure at a point and the slope of the body 
surface at that point. In a similar manner, the interference, drag coef- 
ficient for the case when each body lies within the disturbed flow field 
of the other can be calculated. 

As a preliminary step, therefore, it is necessary to find an expres- 
sion for the pressure around each body for the whole region behind the 
Mach wave from the nose of the body. 
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METHOD OF CALCULATION 

Tne Pressure Field Surrounding A Single Body 


Under the assumptions of linearized theory, the shape. of a slender 
body of revolution is described by the distribution along an axis of 
sources and sinks which satisfy the potential equation and the boundary 
conditions of uniform flow at infinity. Under these assumptions the 
source strength is given by 


2*f’(x) = V- 


d^ 

dx 5 


where S is the expression for the area of a section. 


The pressure coefficient, in this theory, is found from the relation 

( 2 ) 


c p . _ - 
* V dx 


where, as given in reference 2, d<p/dx is given by 

nose J (x-£) 2 ~ 


<yp = _ r* 
“ J 


( 3 ) 


Here 


f'(l) - 


m2 - 2(S-b) 


2 1 


-J m 2 - ( 5 -b ) 2 


is the expression for the source strength in terms of the following 
coordinates of integration: 


B 2 = M 2 - 1 


P P P 

r = y + z 


When f * ( |) is replaced in equation ( 3 ) by its equivalent, and the 
substitution is made for d9/dx in equation (2), the pressure coeffi- 
cient is found to be determined by the integral 


= - 3C / 


[m 2 - 2( |-b) 2 ]d| 


L J Em 2 - ( l-b) 2 ] [(x-i ) 2 - B 2 r 2 ] 


(k) 
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This integral is elliptic, and the limits of integration and the method 
of solution depend on the following tnree regions (see fig. 2): 

Region I is bounded by the Mach aftercone from the nose and 
the Mach forecone from the tail. The limits of integra- 
tion are b-m and x-Br. 

Region II is bounded by the Mach forecone from the tail and 
the Mach aftercone from the tail. The limits of integration 
again are b-m and x-3r. 

Region III is bounded by the Mach aftercone from the tail and 
infinity. The limits of integration are b-m and b+m. 

With these limits the elliptic integrals for the different regions 
are complete and the solutions are: 1 


Region I 



2 J (m+3r) 2 - (x-b) 2 
[(m+3r) 2 - (x-b) 2 ] E 0 (kl) + 



(5) 


where 



i+3r+x-b 

2m 


and the change of variable 


J (m-3r-b+x)(m+b-|) 


y (m+3r+b-x) ( |«n-b ) 
(m-3r-b+x)(m+b-|) 


(m+3r+b-x) ( l-nn-b) 


transforms equation (4) into normal elliptic form. 


^Solution of these integrals was accomplished with the aid of a table of 
elliptic integrals compiled by Mr. Paul F. Byrd of the 
Ames Aeronautical Laboratory. 
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Region II 


= . |m(m+20r+2b-2x)F o (k II )-4m0rE o (k II )+U(x-b)ymSr A 3 (k II ,£ II )J- 

2 J m0r (6) 


where 


kll =. / (^) £ - 5n „ sin -1 /-— 

</ l4m0r v x- 


2m 


•b+Sr+m 


and the change of variable 


u = sn 1 / 

y (*-i 


23 r(l-b+m) 


-b - 0r+m ) ( x+ 0r- 1 ) 
transforms equation (10 into normal elliptic form. 


Region III 


C p = 


+3*c 


(x-b-3r) 2 F 0 (kux) + [(x-b) 2 - (m-3r) 2 ] 


*/(x-b) 2 - (m-3r) 2 
[Eo^IIl) 3 “ 2(x-b) J (x-b) 2 - (m-3r ) 2 A 0 (k III , ^ m ) j- 


where 


k TTT = / > Bm .1; U = sln-i /EEid 
111 y (x-b) - (m-0r) 2 J x+Sr+m- 1 

and the change of variable 


_i / (x+0r-m 
u = sn 1 / 1 - 

V 2m(x 


-m-b)(i-b-m) 


2m(x+3r-| ) 

transforms equation ( k) into normal elliptic form. 


(7) 
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The functions Fo> Eo> and A 0 are tabulated in reference K, where 
they are also defined as: 


FoOO = -|K(k) 

E 0 (k) = -f-E(k) 


Ao(k,C) = E 0 (k)F(k’ ,C) + F 0 (k)E(k»,£) - F 0 (k)F(k»,S) 


where 

K(k) 

E(k) 


complete elliptic integral of the first kind of modulus k 
complete elliptic integral of the second kind of modulus k 


F(k* Cll 

E(k* * £)j i ncom Pl e ^ e ellipti c inte grals of the first and second kind of 
modulus k* (= »J 1-k 2 ) , and argument £, respectively 


Thus, according to linearized theory, equations (5), (6), and ( 7 ) 
completely determine the pressure field in the region behind the Mach 
wave from the nose of a body of revolution of the prescribed shape. An 
isometric sketch of the pressure coefficient at fixed values of the 
radial distance is shown in figure 3* 


Interference Drag Coefficient 


An inspection of figure 3 shows that the gradient of the curve of 
the pressure coefficient in the stream direction changes from negative 
to positive behind the center of the body. This region of positive 
pressure gradients is a zone of favorable buoyancy in the pressure field 
around the body. Therefore, a small body placed anywhere within this 
zone should have a negative interference drag (i.e., a thrust) due to 
this pressure field which should cause a decrease in the combined drag. 


If the pressure coefficient surrounding a body at a given lateral 
distance r = ro be called Cp^, the interference drag coefficient may 
then be calculated from the expression 


Cdj - 


rtRmax ‘ 


/ 

'-'nr 


tail 


'Pt 


nose 


_ RdR 

2n dx 

dx 


(8) 
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where Rmax is the radius of the frontal area of the reference body. 
Substitution of the respective values for the quantities in equation ( 8 ) 
yields 


Actually, the integral in equation (8) or ( 9 ) is made up of one or more 
integrals depending on the number of regions in which the small body is 
located. The limits of integration depend on the parameters b, m, r, 
and 3 which determine the regions of integration. This integral must 
be computed numerically. 

For a clearer exposition of the mathematical computations in the 
cases where interference exists, it is necessary to distinguish between 
the bodies. In the present case, if the' central body situated on the 
x axis be denoted by the subscript 1 and the other body by the sub- 
script 2 , then, depending on which body is being acted upon, the inter- 
ference drag coefficient is 


v Equations (10) and (ll) are, to a first approximation, the drag 
coefficients due to the effect of body 1 on 2 and body 2 on 1 , respec- 
tively. The. integration extends from the nose to the tail of the body 
which is in the disturbed flow field. Similarly, for the calculation, 
to a first approximation, of the drag coefficients due to the inter- 
action of the bodies, formulas ( 10 ) and (ll) are both used. 

The preceding equations permit the calculation of the interference 
drag coefficient for any arrangement of bodies. In particular, equa- 
tions ( 10 ) and (ll) apply to the special case when two identical bodies 
are symmetrically placed with respect to a central body so that each 
body may be in the flow field of the others. The formulas may be applied 
to any configuration, whatever the interference pattern. 



Cp t ^ - 3 *c (x-b) J m 2 - (x-b) 2 j> dx (9) 



(10) 


or 
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Numerical Calculations 


As an exploratory investigation, the simple case of a small body 
having one-half the length of the large body was chosen. This combina- 
tion is derived by assigning the following values to the parameters in 
equation (l): 


It is to be noted that each body has the fineness ratio 10 but that one 
is twice as long as the other. , , 

The pressure coefficient depends on 3 and r 0 . 'The dependence is 
on 3r 0 , rather than on either parameter individually, and the calcula- 
tions were carried out for 0r o = 0.5m 2 , 1.0m 2 , and 2.0m 2 ; where^ as 
shown previously m 2 = 1. A change in 3 and in ro which keeps 3i*o 
invariant will not change the calculations. 

The drae coefficient, with resnect to the frontal area, of the 

in 


formula ( 9 ). For each value of 3r 0 , therefore, the remaining parame- 
ter b, in this case the longitudinal distance between the centers of 
the bodies, is given values between -3m 2 and +4m 2 . 

The interference drag coefficients based on frontal area of the 


the reversibility of drag (reference 5) holds true, the curves would be 


nondimens ional longitudinal distance between the centers of the bodies). 
The symmetry in the figure is a measure of the accuracy of the numerical 
calculations. From the figure, it is evident that minimum drag occurs 


when the small body is close to the large one and with 

its center just forward of the tail wave from the large body. 


c 


m 


b 


C]_ = 0.005 
c 2 = .01 


mi = 2 b 1 = 0 

m 2 — 1 b 2 = b 


The equations which describe the bodies are 


R x 2 = 0.005 (4-x 2 ) 3/2 
R 2 2 = .01 [l - (x-b) 2 ] 3/2 


I 



small body are sketched for different and — in figure 4. Since 

I H2 ^2 


expected to be symmetrical about the line — = 0 (where ~ is the 

IH2 Bljg 
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It may be observed that, while the interference drag varies markedly 
with fore and aft position, for the range of pr Q considered the maxi- 
mum favorable interference that can be obtained generally decreases with 
increasing separation between the bodies. It is interesting to note 
that drag minimums occur whenever the center of the satellite body is 
just forward of the tail Mach -wave of the parent body. Also of interest 
is the fact that, in the case when both bodies had their centers on the 
vertical axis, the lowest drag occurred when the lateral separation 
parameter (pro) was equal to the length of the small body and higher 
drags resulted as the bodies were brought closer together. 

In figure 5# there are plotted the values of the total. drag coeffi- 
cient, based on the total frontal area of a three-body combination with 
bilateral symmetry, against the longitudinal distance between the centers 
of the parent and satellite bodies. Since interaction between all three 
bodies occurred only for the lowest value of Pr© (= 0 . 5 ^ 2 ) where it was 
found to be negligible, the interference drag coefficients previously 
calculated for the two-body configuration could be used directly to 
determine the total drag of the three-body configuration. As a conse- 
quence, the variations in drag coefficient are similar to those of 
figure 1+. Again the lowest drag occurred at pro = 0.5 m 2 an d was 
approximately 35 percent less than the drag of the three bodies without 
interference. 

Because of the unusual shapes of the curves in figure k, it was 
considered advisable to investigate the drag interaction of a combina- 
tion of bodies of different shape. The slender pointed body derived by 
Jones and Margolis (reference 6) was selected for this purpose since for 
the same fineness ratio the drag coefficient is comparable to that of 
the Sears-Haack body. The interference drag coefficients of a combina- 
tion of two such bodies with mi/mg = 2 were calculated for b = 0, and 
different values of Pr D . Reasonable agreement was obtained with the 
results of figure 4 for pr Q = 2.C*n 2 and 3r 0 = 1.0m 2 . In the case when 
Pr 0 = 0.*m 2 , since the bodies are close together, there is a discrepancy 
which may be due to the differences in body geometry. , 


CONCLUDING REMARKS 


It is found that the combined wave drag of a combination of bodies 
of revolution can be decreased if an arrangement is chosen which takes 
advantage of a favorable pressure zone which exists behind the center of 
each body. 

In the cases discussed, where the ratio of the lengths of the 
bodies was 2 to 1 but the fineness ratios were equal, numerical calcula- 
tions showed that the maximum favorable interference occurred when the 
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center of the small body was just forward of the stern Mach wave from 
the large body and the bodies were close together. For the range of 
Pr Q considered the magnitude of the favorable interference generally 
decreased with increasing separation between the bodies. 

In the case of a bilaterally symmetrical arrangement of three bodies 
with a lateral separation equal to one -quarter of the length of the small 
body the total wave drag was found to be 35 percent less than the com- 
bined wave drag of the three bodies. 


National Advisory Committee for Aeronautics 
Ames Aeronautical Laboratory 

Moffett Field, Calif., Sept. 20, 1951 
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Figure /.- Coordinate system. 
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Figure 2.- Regions of integration. 
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Figure J.- Isometric sketch of the pressure distribution in the field surrounding a single body. 
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Figure 4.- Interference drag coefficient of configuration based on the area of the small body. 
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Figure 5.- Total drag coefficient of three-body combination with bilateral symmetry . 
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